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13.  Abstract 

A  new  theoretical  model  for  the  development  of  an  estuary  bed  profile  under  the  influence  of 
river  flow  is  presented.  The  model  is  based  on  original  dependences  describing  the  distribution  of  the  free¬ 
water  surface  slope  both  along  and  across  the  mouth  reach  that  were  developed  by  an  analysis  of  experi¬ 
mental  data.  A  fundamentally  new  solution  is  given,  which  accounts  for  the  turbulent  interactions  between 
flow  jets  in  the  coastal  area  where  the  river  flow  actively  interacts  with  the  surrounding  water. 

Two  variants  of  the  theoretical  model  were  formulated,  a  one-dimensional  and  a  two- 
dimensional  model  (a  flow  plan).  Numerical  solution  of  the  one-dimensional  model  allowed  prediction  of 
bottom  deformations  with  time  under  the  influence  of  key  hydrological  factors  at  different  slopes  of 
coastal  bottom  relief.  Numerical  solution  of  the  flow  plan  model  provided  the  longitudinal  and  transverse 
distributions  of  flow  velocities  and  bottom  relief  deformations. 

An  analysis  of  the  obtained  information  afforded  novel  dependences  describing  the  rates  of  bot¬ 
tom  relief  formation  and  bar  crest  growth,  as  well  as  the  final  height  of  the  bar  crest  and  the  rate  of  its 
downstream  migration.  The  results  of  calculations  are  in  perfect  agreement  with  the  experimental  data  and 
field  measurements.  The  results  of  this  work  will  help  choosing  the  right  location  and  depth  of  navigable 
channels  in  estuaries. 

14.  Mouth  bar,  Estuary,  Navigation  draft,  Numerical  simulation 

15. 19  Pages 

16. 

17.  Unclassified 


18  Unclassified 


19.  Unclassified 


20. 


Table  of  Contents 


Technical  Objective  2 

Background  2 

The  theory  of  the  problem  2 

The  behaviour  of  the  free-water  surface  slope  3 

Flow  velocity  within  the  estuary  reach  4 

Flow  velocity  within  the  coastal  area  5 

Bottom  relief  deformations  9 

Method  of  calculations  9 

The  results  of  calculations  1 1 

Distributions  of  flow  velocities  1 1 

Bottom  relief  deformations  1 1 

Conclusions  12 

References  1 3 

Appendix  1.  Illustrations  1-1 


1 


Technical  Objective 

The  objective  of  this  work  was  to  create  a  universal  theoretical  description  of  the  develop¬ 
ment  and  deformations  of  mouth  bars  caused  by  a  combination  of  river  and  sea  hydrological  factors. 
The  developed  model  assumes  the  absence  of  considerable  tidal  and  sea-wave  influence  on  the  mouth 
bar  relief:  Development  of  such  a  theory  should  be  considered  as  a  substantial  contribution  to  hydrol¬ 
ogy. 

The  proposed  numerical  model  allows  prediction  of  1)  mouth  bar  deformations  due  to  sea¬ 
sonal  changes  in  the  hydrological  regimes  of  both  river  and  sea  and  2)  obstruction  of  the  navigable 
channel  under  the  influence  of  active  factors  in  various  combinations.  Successful  solution  of  this  prob¬ 
lem  opens  new  horizons  in  the  prediction  of  mouth  bar  deformations. 

Background 

The  cyclic  development  of  delta  branching  is  extensively  studied  by  drilling  and  documented 
in  archive  sources  [11  The  redistribution  of  flow  in  delta  branches  has  an  adequate  theoretical  descrip¬ 
tion  [2,3].  Unlike  these  processes,  the  morphological  changes  in  bar  relief  that  occur  under  the  com¬ 
bined  influence  of  river  and  sea  hydrological  factors  have  no  analytical  solution  as  yet.  Only  a  few 
rough  approximations  such  as  the  so-cailed  "mouth  elongation"  can  be  found  in  the  literature  [4].  Our 
extensive  experimental  studies  [5,6,7],  examination  of  actual  data  on  numerous  river  estuaries  (in  par¬ 
ticular,  those  of  rivers  flowing  into  the  marginal  seas  of  the  Arctic  Ocean  [7;&]Y  and  theoretical  calcu¬ 
lations  [9,10,1 1]  constituted  the  elements  of  a  universal  theory  of  mouth  bar  development  [7,12,13]. 
Fragmentary  realisations  of  this  theoretical  model  were  successful  [14,15,16]. 

The  theory  of  the  problem 

Navigational  traffic  in  river  estuaries  requires  practical  knowledge  about  the  intensity  of 
sedimentation  in  the  estuary,  mouth  bar  development,  navigational  depth  at  the  bar  crest,  and  the  rate 
of  migration  of  bar  elements.  In  the  case  being  considered,  the  navigation-affecting  parameters  of  a 
mouth  bar  totally  depend  on  river  factors.  The  mouth  bar  development  mainly  occurs  at  high-water 
periods  (floods),  when  the  bar  crest  reaches  its  maximum  height  and  the  bar  itself  is  maximum  shifted 
seawards.  This  time  is  least  favourable  for  navigation.  Sea  factors,  especially  wind-induced  waves, 
usually  destroy  the  bar  relief  formed  by  river  flow,  restructure  it,  and  lower  the  bar  crest.  Such  defor¬ 
mations  often  increase  navigational  depths.  However,  the  effect  of  sea  factors  on  the  stability  and  pro¬ 
file  of  a  navigable  channel  is  ambiguous.  The  river-  and  sea-related  processes,  which  are  often  coun¬ 
teracting,  will  be  discussed  elsewhere. 

Variations  of  the  bar  crest  height  with  time,  the  rate  of  its  vertical  deformation  and  longitudi¬ 
nal  migration,  as  well  as  other  parameters  that  depend  on  annual  and  perennial  changes  in  river  runoff, 
ffee-water  surface  slope,  and  sediment  discharge  at  the  upper  boundary  of  the  estuary  reach  will  sub¬ 
sequently  be  referred  to  as  the  "mouth  bar  regime”. 

The  estuary  area  is  usually  subdivided  in  two  regions  [4]:  the  estuary  reach  where  river  fac¬ 
tors  dominate  and  the  coastal  area  with  distinctly  sea-type  hydrological  regime.  The  upper  boundary  of 
the  estuary  reach  normally  coincides  with  either  the  delta  head  or  the  section  line  at  which  no  positive 
set-up  can  be  observed.  The  lower  boundary  is  a  virtual  line  enveloping  either  the  sea  coast  or  delta 
islets  reached  by  the  sea.  Section  line  at  which  the  influence  of  the  river  flow  disappears  separates  the 
coastal  area  from  the  open  sea. 

The  theoretical  model  uses  the  system  of  orthogonal  curvilinear  co-ordinates  shown  in  Fig.  1. 
The  origin  of  co-ordinates  is  in  the  centre  of  the  estuary  section  line,  The  longitudinal  co-ordinate  /  is 
directed  downstream.  The  upper  and  lower  boundaries  of  the  estuary  area  (as  defined  above)  were  as¬ 
sumed  to  be  at  —  oo  and  4-  so,  respectively.  In  the  one-dimensional  model,  flow  velocities  and  bot¬ 
tom  deformations  are  either  considered  along  the  stream  axis  (the  isoline  of  maximum  depths)  or  aver¬ 
aged  across  the  flow,  In  the  two-dimensional  model,  the  river  flow  is  divided  into  separate  flow  jets, 
This  approach  known  as  the  construction  of  a  flow  plan  simplifies  both  the  set-up  and  solution  of  the 
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problem  compared  to  the  use  of  orthogonal  co-ordinates.  The  ordinate  axis  is  denoted  b  axis.  Flow 

velocities  are  averaged  over  depth. 


The  behaviour  of  the  free-water  surface  slope 

It  is  known  that  the  free-water  surface  slope  of  a  river  depends  on  such  parameters  as  bed 
profile,  bed  roughness,  and  flow  rate  and  is  mainly  determined  by  the  surface  drop  between  the  upper 
boundary  of  the  estuary  reach  and  the  reservoir  level.  As  water  level  at  the  upper  boundary  greatly 
varies,  whereas  the  sea  level  is  constant,  the  maximum  surface  slope  should  correspond  to  catastrophic 
floods,  and  the  minimum  slope,  to  low-water  periods. 

River  flow  above  the  estuary  is  usually  considered  to  be  quasi-uniform,  with  a  constant  sur¬ 
face  slope  along  the  river  bed.  Within  the  estuary  reach,  river  flow  is  always  non-uniform.  Here  the 
surface  slope  is  increasing  downstream  at  high-water  periods  and  reaches  its  maximum  at  the  estuary 
section  line  (Fig.  1).  Location  of  the  maximum  slope  depends  on  flood  intensity.  The  higher  the  flood, 
the  more  the  section  line  of  maximum  surface  slope  is  shifted  seawards.  Within  the  coastal  area,  the 
surface  slope  is  first  sharply  decreased  (as  the  river  flow  is  no  longer  confined  between  the  banks)  and 
then  gradually  disappears.  At  low-water  periods,  a  continuous  decrease  of  the  surface  slope  down  the 
river  is  observed  virtually  everywhere  along  the  estuary  reach.  Within  the  coastal  area,  the  surface 
slope  is  only  observable  at  the  bar  crest  during  negative  wind-induced  set-ups. 

A  statistical  analysis  of  experimental  results  and  field  measurements  provided  an  analytical 
expression  of  the  surface  slope  behaviour.  The  slope  depends  on  1)  the  surface  slope  above  the  estuary 

area  /-  00  ,  2)  the  positive  difference  AZj^  between  the  surface  level  at  the  estuary  section  line 

that  would  be  observed  in  the  case  of  quasi-uniform  water  flow  within  the  estuary  reach  and  the  reser¬ 
voir  level  (Fig.  1),  and  3)  the  morphological  index  kj  of  the  river  bed: 

/  =  /_„  +  /_„«/ ekTC'kl  at  T  <0,  (1) 

I  =  I0e~kTclr  MT > 0,  (2) 

where  /q  +  )/_oq  is  the  surface  slope  at  the  estuary  section  line, 


Uj  is  the  slope  amplitude, 

Q  and  Q  are  the  parameters  of  the  estuary  reach  and  coastal  area,  respectively,  that 

determine  the  degree  of  the  surface  slope  transformation; 

kf  =  /  //_oo  is  the  morphological  parameter  of  the  river  bed 

where  Xf  is  the  hydraulic  friction  coefficient, 

5_oo  and  //_  qo  are  the  average  width  and  depth  of  the  bed,  respectively,  above  the 
estuary  area, 

/  =  /  /  B_o0  is  the  relative  distance  from  the  estuary  section  line. 

An  expression  for  the  longitudinal  profile  of  the  free-water  surface  is  found  by  direct  integra¬ 
tion  of  (1)  and  (2): 


%  =  „  T  <  0, 


7w  =  1  B  0  +  al)  e~kTC[i  r 

JW  1—COJJ—CD  J  ^  at  1 


kTCl 


>0 


(3) 
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where  Zur  and  Zur  ^  are  the  free-water  levels  (measured  with  respect  to  the  constant  level  of 

calm  sea)  in  the  river  bed  and  at  the  upper  boundary  of  the  estuary  reach  (the  initial  section), 
respectively. 


The  slope  amplitude  Uj  can  be  determined  using  the  limiting  condition  of  smoothness  (the 

absence  of  discontinuities)  imposed  on  the  curve  that  represents  the  free-water  surface  at  the  estuary 
section  line: 


dZyy 


dl 


_  dZyy 


L»-G 


dl 


or  I  i 


/->+ 0 


:/->-0  “  I/-M-0' 


Tills  condition  together  vvitli  ( 1  ^  snd  (2^  gives  the  following  cuprcssion  for  the  slope  amplitude! 
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whence  AZj  j/j  ,  which  is  the  sum  of  free-water  level  drops  within  the  estuary  reach  and  coastal  area 
(Fig.  1 ),  is  determined: 


/_*>£_ 


A  7  _[  g/  l  +  aj 

w°  lq  q  j  kr 


QO^— 00 


r^o  theoretical  justification  of  the  C fa  and  values  exists  as  yet.  In  the  present  model, 

these  parameters  are  free  variables  determined  by  comparing  the  experimental  data  with  the  results  of 
numerical  modelling.  Either  parameter  may  be  greater  or  less  than  unity.  An  analysis  of  experimental 
data  and  field  measurements  showed  that  the  curve  representing  the  free-water  surface  is  smoothed 
during  the  development  of  a  mouth  scour  and  formation  of  the  bar  crest  with  both  the  discharge  and 
level  of  water  in  the  initial  section  being  constant.  This  smoothing  can  be  accounted  for  by  decreasing 

the  Cfc  parameters.  A  simple  empirical  formula  was  discovered: 

q-q-q-i/tfmax. 

where  //nax  ~  ^Ana\.  ^  co  > 

is  the  maximum  depth  of  mouth  scour,  and 
/L  oo  ls  the  river  depth  above  the  initial  section. 

Assuming  that  ^~'jc  ^  at  t  —  0,  formula  (4)  can  be  transformed  to  express  the 

slope  amplitude  as  a  function  of  time: 

«/(0  =  0.5[(l  +  0.5a/(0))q  - 1], 

It  should  be  noted  that  the  slope  amplitude  Ctj  and  parameters  CL  and  CX  may  differ  for 

each  of  the  jets  in  the  flow  plan  model,  in  general,  they  are  gradually  decreasing  upon  transition  to¬ 
wards  peripheral  jets. 


Flow  velocity  within  the  estuary  reach 

Water  flow  within  the  estuary  reach  is  both  non-uniform  and  unsteady.  However,  the  per¬ 
formed  tests  and  comparisons  (e.g.,  see  [17])  show  that  the  term  accounting  for  the  local  acceleration 
in  the  differential  equation  of  motion  is  negligible  compared  to  other  terms  even  during  build-ups  and 
recessions  of  natural  floods.  It  was  shown  that  the  turbulent  tangential  stress  can  also  be  neglected  in 
this  case  [10].  Therefore,  the  initial  equation  can  be  written  as  follows: 
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d\Vl\  XT  V* 


/  =  —  —  +  ^~  —  at/<0.  (5) 

dl{2g)  H  2g 

Here,  V  and  H  are  the  flow  velocity  and  depth  for  a  given  cross-section.  The  solution  of  ($)  can  be 
written  in  this  form  [5]; 

/  ^ 


2  f  ~^7 

V2  =2g  f  Ie  * 


For  an  open  bed,  the  hydraulic  friction  coefficient  Xf  can  be  expressed  via  the  Chezy  friction  factor 

2 

C  as  Xf  —  2g  /  C  .  Using  the  Manning  formula  for  the  Chezy  factor  and  assuming  a  constant 

coefficient  of  bed  roughness  along  the  flow,  we  obtain  the  following  equation,  in  which  the  hydraulic 
friction  coefficient  is  expressed  via  its  value  in  the  initial  section  (viz.,  at  the  upper  boundary  of  the 
estuary  reach): 

J  f  dl  =  XT_xH’Xl  j 

-00  -  00  ** 

where  X j1  ^  is  the  hydraulic  friction  coefficient  for  the  initial  section  and  TH\  —  4  /  3  . 

Substitution  of  this  expression  and  that  of  the  surface  slope  (1)  in  (6)  followed  by  reduction 
to  a  non-dimensional  form  gives 

U  =  A'l J  'je-^'^dX  +  a,  J  (7) 


where  U  — 


vL  '  ’  lHm' 


—  ,and/?  =  — 


As  water  flow  in  the  initial  section  is  assumed  to  be  quasi-uniform, 

_  ^gl  —  CO B—cc  __  | 

vl^T 

The  obtained  expression  (7)  is  valid  both  in  the  one-dimensional  and  two-dimensional  (flow 
plan)  models. 

Let  us  consider  the  initial  stage  of  bottom  relief  formation,  when  the  depth  is  constant  along 
the  river  bed.  In  this  case,  (7)  takes  this  simple  form: 

U  =  \  +  -^—ekTc'kT  (8) 

i+Q 


Flow  velocity  within  the  coastal  area 

Within  the  coastal  area,  the  river  flow  spreads  unconfmed  by  the  banks,  which  results  in  ac¬ 
tive  turbulent  interactions  between  the  river  flow  and  the  surrounding  water  mass.  This  process,  as 
well  as  bottom  friction,  efficiently  quenches  the  flow  energy  and  decreases  its  velocity.  In  this  case, 
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equation  (5)  should  be  complemented  with  a  term  accounting  for  the  tangential  stress  on  the  lateral 
surfaces  of  flow  jets  [10]. 

Methods  of  calculating  breakdown  flows  applicable  to  fiver  flows  as  well  were  developed  by 
numerous  researchers  [18,19,20,21,22].  Vast  experimental  data  on  the  expansion  of  breakdown  flows 
can  be  found  in  [19,231  A  solution  of  the  problem  accounting  for  the  bottom  friction  and  turbulent 
tangential  stress  in  the  ease  of  both  ffee-water  surface  slope  and  depth  varying  along  the  river  was  re¬ 
ported  in  our  earlier  works  [9,10,24]. 

Given  a  preferred  flow  direetion5  so  that  the  velocity  vector  is  practically  collinear  with  the 
longitudinal  co-ordinate,  the  simplest  expression  for  the  turbulent  tangential  stress  is 

w 


where  the  turbulent  viscosity  coefficient  Vf  is  a  scalar. 

The  up-to-date  models  of  turbulence,  e.g.  the  differential  K  ~  8  model,  seem  to  be  the  best 
founded  from  the  theoretical  viewpoint.  However,  application  of  such  models  to  the  discussed  prob¬ 
lem  would  have  considerably  complicated  the  analytical  expressions  for  stream  velocities  with  only  a 
small  gain  in  accuracy.  Therefore,  we  favoured  the  hypothesis  of  V,M>  Makkaveev  [25]  that  the  tur¬ 
bulent  viscosity  coefficient  is  proportional  to  averaged  velocity,  viz. 

Vt-AV, 

where  A  is  a  function  with  linear  dimensionality. 


To  deduce  an  equation  for  the  velocity  distribution  in  a  free  turbulent  jet,  LM.  Konovalov 
[26,27]  treated  the  diameter  (width)  of  the  expanding  jet  as  a  characteristic  dimension  and  expressed  it 

as  a  linear  function  of  the  longitudinal  co-ordinate.  The  resulting  expression  for  Vf  was 

vt  =  2a1  IV  (9) 

where  a  is  the  turbulent  exchange  coefficient. 


A  similar  dependence  was  earlier  discovered  by  Reichardt  [28]  who  analysed  the  experimen¬ 
tal  data  on  the  velocity  distribution  in  a  turbulent  jet.  The  phenomenological  theory  of  Reichardt  at¬ 
tracted  serious  criticism.  Nevertheless,  it  is  in  a  good  agreement  with  the  experimental  data  and  has 
undeniable  practical  advantages  [29] .  Assuming  the  Reichardt— Konovalov  definition  of  the  turbulent 
viscosity  coefficient,  we  arrive  at  the  Following  differential  equation  of  flow  motion  that  includes  all 
the  key  factors  of  turbulence; 
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Substituting  V~  —  ue  ^  J  —  —  — —£  ^  an(j  qp  —  u!%g 
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obtains 
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4-  Z .  one 


(10) 


As  river  flow  practically  expands  in  the  coastal  area  as  in  an  unbounded  half-space,  the  initial  condi¬ 
tions  are 

<4=0  =  =  v0  /2& + *0 ,  if  -  Bo  1 2  </><  Bo  [2 . 
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,  f  b>B0/2 

l/=0  \b<-B0/2’ 

where  Vq  and  Zq  are  the  flow  velocity  and  the  Z  value,  respectively,  in  the  mouth  section. 


Following  is  the  solution  of  (10): 


2o14tc 


f  *(«)< 


\2  al) 


- Bn/2 


Flow  velocity  can  be  written  as 


V2  =  V$IW+2 gz0(ZW - 1)  +  2g\ Ie^F{l)dl  e  V{1\ 


where  E1?  =  —  erf{- 1/2 +  ^ 

2  V  2 al/B0 


1/2 

2aliB0 


2  xr 


evf  (V)  —  —j=  J  e  t  dt  is  the  error  function, 


Bq  and  Hq  are  the  bed  width  and  depth,  respectively,  in  the  mouth  section, 

X is  the  hydraulic  friction  coefficient  in  the  mouth  section. 

Changing  over  to  dimensionless  values  and  taking  account  of  the  expression  (2)  for  the  sur¬ 
face  slope,  one  can  write 

u  =  + (Uoo  -  Uofev  - 1 )Ykr^'  -°)  + 

+  (1- +  a,  )kT \fkr  K'  ’*)+C**W 


where  Uq  =  — y — ,  Uqo  =  -y®- 

VL.O  V-m. 0 


Vqq  and  are  the  axial  flow  velocities  in  the  mouth  section  and  at  the  upper  bound¬ 

ary  of  the  estuary  reach,  respectively. 

The  solution  (12)  is  valid  at  any  point  along  the  river  flow  where  the  confining  banks  do  not 
substantially  affect  the  flow  motion.  In  the  particular  case  of  the  one-dimensional  model,  flow  velocity 
is  described  by  the  following  expression: 
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U  =  U( 


J  6a  \ 

whereSf- . 

At  the  initial  stage  of  bottom  relief  formation,  the  depth  is  everywhere  constant,  therefore 

t/s  t/oE'F  +  {1 

_  '"'fc  _ 

There  is  an  alternative  non-dimensional  representation  of  (1 1)  for  the  case  of  an  axial  jet  ; 

4  =  +  AnAM-^^cR, 


where  A”  —  . 

vikT 

Let  us  consider  the  physical  meaning  of  parameter  ^4”  .  As  was  shown  above,  ^4*  =  1  in 
the  initial  section  regardless  of  water  level  In  the  mouth  section,  Arf  can  be  either  greater  or  less 
than  unity  5  le:  A”  >  1  during  floods  and  A”  <  1  during  low-water  periods.-  It  follows  from  (3) 

iKni  iUa  «\/\o>+ixja  rlfPPaf  ka+xtjiiar*  iUa  fVaa.ttjo+iir  tmrfnrta  1  ax  ml  in  ika  mrviiiU  p  n  nt  i  i-v«  nnrl  ika  ciakla  aaa 

tual  tiiv  puom  vw  \1iaIvi  uwi  vvwui  Uiv  jli  v>C-  vv a.tCi  oui  idtv  ibvvi  In  liiv  mv^utii  oCuliUil  cuiva  tuv  3t«.uiv  oCa 


Wl  % tC'/c  krC'k 

Substitution  of  this  formula  in  the  above  expression  for  A”  gives: 


IqB-oo  _  (l  +  aj)l_o0B_ 


4"  =  T2  VCT 

(Vol2g) 

This  is  a  ratio  between  the  excess  specific  potential  energy  of  the  flow  in  the  mouth  section  with  re¬ 
spect  to  the  reservoir  level  and  the  specific  kinetic  energy  of  the  flow  in  the  same  section  multiplied  by 

Q. 

Assuming  Cj£  =  I  >  one  arrives  at  a  simple  relationship  between  the  distribution  of  velocities 

along  the  coastal  area  and  the  quotient  A Zffiy  /  y^Q  /  Thus,  when  the  potential  energy  of 

the  flow  jet  is  greater  than  its  kinetic  energy,  there  occurs  a  partial  transition  of  the  former  into  the 
latter.  Accordingly,  the  flow  is  accelerated,  and  its  velocity  reaches  a  maximum  at  some  distance 
downstream  from  the  mouth  section,  Then,  the  flow  velocity  gradually  decreases  due  to  turbulent  and 
frictional  dissipation  of  energy.  If  the  potential  and  kinetic  energies  are  equal,  the  flow  velocity  slowly 

dV 

decreases  along  the  longitudinal  co-ordinate,  so  that  — —  —  U .  If  the  potential  energy  is  less 

dl  /- o 

than  the  kinetic  energy,  the  flow  jet  velocity  abruptly  falls  as  soon  as  the  flow  reaches  the  coastal  area. 
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Bottom  relief  deformations 

Bottom  deformations  of  the  estuary  reach  and  mouth  bar  depend  on  the  actual  discharge  of 
bed  load,  which  is  determined  by  the  hydraulic  characteristics  of  the  flow  in  any  cross-section  of  the 
river  bed.  Our  investigations  showed  the  distribution  of  sediments  along  the  river  to  be  a  function  of 
flow  velocity.  The  streamwise  migration  of  sediments  is  practically  unaffected  by  depth  changes.  This 
together  with  the  fact  that  the  water  discharge  is  constant  along  the  flow  trajectory  gives  an  expression 
for  the  rate  of  vertical  deformations  of  the  river  bed  and  mouth  bar: 

dH  qs_„  (  U  dH  3  af\ 
b  dt  (i - s)B_^  \h a  2 a) 

where  qs  is  the  sediment  discharge  in  the  initial  section  (volume  units  per  unit  of  bed  width) 


£  is  the  porosity  coefficient  of  bottom  sediments. 

The  maximum  amount  of  bed  load  that  can  be  displaced  by  river  flow  under  given  hydrologi¬ 
cal  conditions  is  determined  using  the  following  original  formula: 

=  0.0065^1  4dvS— )  ,  (16) 


where  Vn  is  the  non-displacing  velocity  defined  by  V.N.  Goncharov  [30]  as 

Vn  =  \.25jgdlog(8.SH/d), 

V  is  the  kinematic  viscosity,  and 

d  is  the  size  of  sediment  particles. 

The  actual  discharge  of  bed  load  during  different  seasons  of  hydrological  year  is  usually 
smaller  than  the  amount  determined  from  (16).  Therefore  a  correction  factor  should  be  introduced  in 
the  formula  of  the  deformation  rate  to  obtain  correct  quantitative  results.  This  correction  factor  can  be 
estimated  by  comparing  the  actual  and  calculated  data  for  high-water  periods. 

Assuming  that  the  development  of  the  longitudinal  bed  profile  in  the  estuary  reach  and  of  the 
mouth  bar  both  start  when  the  depth  is  still  constant  along  the  estuary  area,  the  initial  rate  of  bottom 
deformations  is 


l-5gj_a 
(1  ~s)B. 


Method  of  calculations 

Since  river  depth  depends  on  the  flow  velocity,  which  in  turn  depends  on  the  occurring  bot¬ 
tom  deformations  and  varies  with  both  water  level  and  runoff,  it  was  necessary  to  divide  the  period  of 
time,  during  which  the  mouth  bar  is  developed,  into  short  time  intervals  such  that  flow  velocity  could 
be  assumed  constant  during  each  interval.  On  every  next  time  interval,  the  flow  velocity  was  recalcu¬ 
lated  to  account  for  the  altered  bottom  relief,  water  level,  and  runoff. 

Calculations  were  performed  using  the  finite  differences  scheme,  which  provided  stable  solu¬ 
tions,  as  well  as  the  required  accuracy. 

In  the  one-dimensional  model,  the  deformation  equation  (15)  was  rewritten  in  this  form: 
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H(t  +  A tj)  -  H(t,l)  -  EU(t,l)  x 

where  A/  is  the  computation  step  aiong  the  longitudinal  co-ordinate  and  A t  is  the  time  step. 

The  dimensionless  quotient  E  is  determined  from  the  initial  parameters  at  the  upper  bound¬ 
ary  of  the  estuary  reach: 

E  _  hsqM  _  Ar 
(\-€)b^h^m  l  a/ 

A  stable  solution  requires  that  the  parameter  E  in  (18)  was  sufficiently  small.  For  the  actu¬ 
ally  observed  flow  velocities  and  depths,  the  condition  E  <  1/4  is  sufficient. 

Calculations  of  bottom  deformations  in  the  two-dimensional  (flow  plan)  model  are  more 
complicated,  Variations  of  jet  width  due  to  changing  depth  and  flow  velocity  requires  to  calculate  a  jet 
depth  taking  into  account  the  depth  of  the  neighbouring  jets. 

The  relative  width  of  each  jet  is 

Bm  =  \/MHmy[U m 

where  Bm  =  Bmj  Bq\  is  the  relative  width  of  an  TH  4h  jet5 
Bm  is  the  actual  width  of  this  jet,  and 
M  is  the  number  of  jets  in  the  plan. 

q  m-\  q 

The  distance  between  the  flow  axis  and  the  axis  of  the  Hi  -th  jet  is  bm  —  -f  ^  Bfc  -h  . 

2  k=2  2 

A  new  depth  of  the  jet  H^{t  +  Af  ^  in  the  same  boundaries  is  determined  from  (18).  A  change  in 
the  jet  width  will  be  then 


A Brn{t  +  At)  = 


V  ^  +  Ai)Vm(t  M &m{i)  ABm^\{t  +  At)] 


e  we  assume  that  H'm+ 1  -  =  1  outside  the  jet  boundaries  and  ABq  =■  0.  New  depths 
2  ' 

r  boundaries  will  be 


jj  ^  \  At)  +  ABm{t  4-  At)Hm+i(t  4-  At) 

Bm(t)+  ABm(t  +  At) 

For  numerical  calculations,  the  integrals  appearing  in  (7)and  (13)  (the  one-dimensional 
model)  or  in  (7)  and  (12)  (the  two-dimensionai  model)  were  replaced  by  the  corresponding  sums,  and 

the  velocity  function  U  was  calculated  using  the  values  of  hydraulic  parameters  II  and  Cfc  corre¬ 
sponding  to  the  end  of  the  preceding  time  interval.  As  there  appears  an  infinite  integral  in  (7),  summa¬ 
tion  started  at  /  ~  —10. —20  The  velocity  and  depth  upstream  from  this  point  were  calculated 
using  (8)  and  (17),  respectively. 
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The  results  of  calculations 

Distributions  of  flow  velocities 

Numerous  experiments  were  performed  to  study  the  velocity  field  in  a  flow  expansion  area. 
However,  most  of  them  either  do  not  conform  to  our  theoretical  model  or  provide  insufficient  data. 
Thus,  data  on  the  local  free-water  surface  slopes  that  are  crucial  for  our  model  are  absent  from  most 
of  these  reports.  Therefore,  we  carried  out  special  experiments  to  establish  the  effect  of  bed  roughness 
on  the  distribution  of  velocities  (averaged  over  depths)  in  the  flow  plan  (the  experiments  used  two 

types  of  bottoms  with  the  hydraulic  coefficients  X]^  =  0.008  and  Xj1^  =  0.014)  and  to  sub¬ 
stantiate  the  values  of  the  turbulence  coefficient  CL  accepted  in  the  model. 

Theoretical  and  experimental  epures  of  velocity  distributions  in  the  flow  plan  are  shown  in 
Fig.  2.  The  experimental  epures  correspond  to  velocities  averaged  over  depths.  The  theoretical  epures 
were  calculated  using  (14)  with  various  values  of  the  coefficient  Q .  The  best  fit  is  observed  at 

Cl  =  0.05  and  Cfc  =  1 .  The  observed  agreement  between  the  calculated  and  experimental  epures 

indicates  that  both  the  accepted  phenomenological  theory  of  turbulence  and  generalisations  of  the  type 
(10)  are  valid. 

At  a  >  0.05,  we  observed  an  increase  in  the  flow  expansion  angle  and  significant  devia¬ 
tions  of  the  calculated  epures  from  the  experiment.  Therefore,  our  attempt  to  account  for  the  bottom 
roughness  by  merely  increasing  the  coefficient  (X  should  be  considered  as  a  failure. 

Bottom  relief  deformations 

We  consider  the  development  of  a  mouth  bar  in  a  coastal  area  with  very  small  bottom  slope 
and  assume  both  the  initial  depth  along  the  estuary  reach  and  runoff  to  be  constant.  These  conditions 
are  readily  reproduced  in  laboratory  experiments,  and  sometimes  they  can  be  observed  in  the  nature 
during  the  formation  of  a  new  estuary.  The  results  obtained  by  solving  the  one-dimensional  model  are 
compared  with  the  experiment  in  Fig.  3.  A  good  agreement  is  seen,  which  indicates  validity  of  the  de¬ 
veloped  model. 

An  analysis  of  the  calculational  results  based  on  the  above  assumptions  most  clearly  shows 
how  the  key  factors  affect  the  mouth  bar  regime.  The  distribution  of  flow  velocities  along  the  estuary 
reach  and  at  the  mouth  bar  are  shown  in  Fig.  4.  The  longitudinal  profiles  of  the  estuary  reach  and 
mouth  bar  at  various  stages  of  bottom  relief  formation  are  shown  in  Fig.  5  with  the  key  parameters 
corresponding  to  the  natural  conditions. 

An  analysis  of  the  results  based  on  the  one-dimensional  model  indicates  that  the  maximum 
free-water  surface  slope  is  located  exactly  at  the  mouth  section  in  the  initial  period  of  the  longitudinal 

bottom  profile  formation.  The  greater  the  surface  slope  amplitude  G]  ,  the  more  is  the  location  of 

maximum  flow  velocities  shifted  seawards.  In  the  longitudinal  bottom  profile,  the  areas  of  mouth 
scouring  and  bar  crest  formation  are  discerned.  The  initial  location  of  the  bar  crest  depends  on  both 

kf  and  G]  and  is  shifted  further  from  the  mouth  section  with  decreasing  the  morphometrical  pa¬ 
rameter  kj . 

Flow  velocity  in  the  neighbourhood  of  the  mouth  section  increases  with  time,  and  the  location 
of  maximum  flow  velocity  is  shifted  downstream.  The  mouth  scour  is  formed  at  Cl]  >  0 .  The  rate  of 

its  formation  increases  with  both  kj  and  CL]  .  The  intensity  of  erosion  decreases  with  time.  The  lo¬ 
cation  of  maximum  depth  migrates  seawards  with  the  velocity  that  is  practically  independent  from  time 
and  is  mainly  determined  by  the  slope  amplitude  G]  . 

The  evaluation  of  mouth  bar  relief  has  two  stages:  1)  the  formation  of  mouth  bar  by  sediment 
deposition  on  the  initially  flat  bottom  of  the  coastal  area  and  2)  the  seaward  migration  of  the  mouth 
bar  thus  formed,  all  its  key  parameters  remaining  the  same.  The  bar  formation  takes  the  time  deter- 
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mined  by  kj  and  CLi .  The  final  relief  is  formed  more  rapidly  with  increasing  kj  and  CLx  parame¬ 
ters. 

Tbs  depth  at  the  crest  of  the  formed  mouth  bar  is  *~Icr  ~~  ^ cr  j oo  =  0.1.  .0.12, 

which  value  weakly  depends  on  kj  and  G]  .  Thus,  the  depth  at  the  bar  erest  is  slightly  decreased 
with  increasing  these  parameters.  The  formed  mouth  bar  migrates  seawards  with  a  constant  velocity, 
which  mainly  depends  on  Qj  and  increases  with  increasing  this  parameter.  An  insignificant  smoothing 
of  the  bar  crest  is  also  observed. 

In  the  course  of  mouth  bar  formation,  the  rate  of  bar  crest  elevation  varies  and  reaches  its 
maximum  when  the  depth  at  the  bar  crest  is  Hcr  =  0.4..0.6 .  The  rate  of  bar  crest  elevation  in¬ 
creases  with  increasing  both  ihe  morphological  parameter  kf  and  slope  amplitude  G  j  .  The  location 
of  the  bar  crest  with  respect  to  the  mouth  section  is  little  changed  during  initial  stages  of  its  formation. 
The  value  of  kj  determines  whether  the  bar  will  migrate  seawards  due  to  sedimentation  on  its  rear 

slope  or  approach  the  mouth  section  due  to  preferred  sedimentation  on  the  front  slope  of  the  bar.  The 
rate  of  seawards  migration  of  the  bar  crest  increases  with  time  and  passes  through  a  maximum  prior  to 
becoming  constant. 

Initially,  the  rear  slope  of  the  mouth  bar  is  very  gentle,  whereas  the  frontal  slope  is  steeper. 
The  steepness  of  both  slopes  increases  with  time,  that  of  the  rear  slope  increasing  more  rapidly.  A 
similar  trend  was  observed  in  our  experiments  [6], 

The  development  of  a  mouth  bar  in  a  coastal  area  sloping  towards  the  sea  is  shown  in  Fig.  6. 
A  comparison  with  Fig.  5  shows  that  the  courses  of  development  and  final  shapes  are  different  for  the 
bars  formed  in  coastal  areas  with  and  without  bottom  slope - 

Calculations  for  the  case  when  floods  of  various  intensity  pass  down  the  river  showed  that 
the  dynamics  of  longitudinal  bar  profile  is  mainly  affected  by  the  level  and  duration  of  floods,  as  well 
as  by  sediment  discharge.  Most  intensive  restructuring  occurs  during  catastrophica!  floods  when  the 
bar  erest  is  noticeably  shifted  seawards.  Normal  floods  (not  to  mention  low-water  periods)  that  follow 
the  catastrophic  ones  practically  do  not  induce  any  significant  bottom  deformations  neither  along  the 
estuary  reach  nor  at  the  mouth  bar. 

The  results  of  flow  plan  calculations  are  shown  in  Fig.  7.  The  initial  bottom  relief  was  a  fiat, 
nearly  horizontal  surface.  In  foil  agreement  with  the  experiment,  maximum  sedimentation  in  the  initial 
period  of  mouth  bar  formation  occurred  in  the  area  of  outer  jets,  where  estuary  sand  bars  were 
formed.  As  the  bottom  is  elevated  with  time,  its  relief  starts  to  affect  the  flow.  The  sand  bars  turn 
around  and  migrate  along  normals  to  the  flow  axis,  which  development  is  in  agreement  with  laboratory 
experiments.  The  bar  crest  is  shifted  seawards.  The  obtained  results  demonstrate  the  similarity  between 
physical  and  numerical  modelling  in  this  case. 

Conclusions 

The  results  of  numerical  modelling  provide  sufficient  information  on  the  mouth  bar  regime, 
which  is  necessary  for  establishing  the  navigational  parameters  of  a  river  estuary. 

The  one-dimensional  model  describes  either  a  non-branched  river  delta  or  the  estuary  of  a 
single  branch.  In  the  latter  case,  the  distributions  of  river  amoff  and  sedimentation  over  ail  delta 
branches,  as  well  as  the  free- water  surface  levels  at  the  heads  of  the  branches  during  each  characteris¬ 
tic  season  of  hydrological  year  should  be  available.  The  described  two-dimensional  model  (flow  plan) 
gives  the  distributions  of  flow  velocities  and  bottom  deformations  both  along  and  across  the  flow.  A 
number  of  characteristics  of  the  mouth  bar  formation  were  .established  that  dep  end  on  the  hydrological 
river  factors  and  the  bottom  relief  of  the  coastal  area. 
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It  should  be  noted  that  the  development  of  a  mouth  bar  in  coastal  areas  with  sufficiently  steep 

bottoms  can  not  be  studied  without  taking  account  of  the  flow  breakdown  from  the  bottom,  which 

phenomenon  deserves  a  separate  theoretical  treatment. 

In  our  opinion,  a  feasible  future  extension  of  this  work  should  be  a  study  of  the  mouth  bar 

formation  and  transformation  under  the  influence  of  wind-induced  waves. 
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Appendix  1.  Illustrations 


Fig.  1 .  A  system  of  co-ordinates  used  in  the  theoretical  model.  The  longitudinal  changes  of  free-water 

surface  slope. 
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Fig.  2.  A  comparison  of  the  calculated  and  experimental  velocity  distributions  along  and  across  the 

coastal  area. 

=  0.014,#  =  0.05  and  Cl  =  1) 


Fig.  3.  A  comparison  of  the  calculated  and  experimental  bottom  deformation. 

( kf  -  0.35,  a  I  =  0.25 , 1  -  /  is  SI  hours,  2  -355  hours) 
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Fig.  4.  The  distribution  of  flow  velocities  along  the  estuary  area. 

(kf  =  0.6,  aj  =  0.5,  E  =  1.0) 


(kf  =  0.6,  Gj  =  0.5,  E  =  1.0) 
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tig.  7. 


Bottom  detormations  along  and  across  the  estuary  area  (the  tlow  plan  model) 

(kT  =  0.6 ,  aj  =  1.0,  E  = 1.0 ,  t  =  1.8 > 


